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ABSTRACT 

We present new and accurate radial velocity (RV) measurements of luminous stars of all ages (old 
horizontal branch, intermediate-age red clump, and young blue plume, as well as red giants of a range 
of ages; 20. 6< V<22) in the Carina dwarf spheroidal galaxy, based on low-resolution spectra collected 
with the FORS2 multi-object slit spectrograph at the VLT. This data set was complemented by RV 
measurements based on medium and high-resolution archive spectra of brighter (V<20.6) Carina 
targets collected with the GIRAFFE multi-object fiber spectrograph at the VLT. The combined 
sample includes more than 21,340 individual spectra of «2,000 stars covering the entire body of the 
galaxy. The mean (<RV>=220.4±0.1 km s^ 1 ) and the dispersion (<r=11.7±0.1 km s" 1 ) of the RV 
distribution of candidate Carina stars (^1,210 objects, 180<i?V~<260 km s -1 , 4er) agree quite well 
with similar measurements available in the literature. To further improve the statistics, the accurate 
RV measurements recently provided by Walker et al. (2007) were also added to the current data set. 
We ended up with a sample of ^1,370 RV measurements of candidate Carina stars that is «75% larger 
than any previous Carina RV sample. We found that the hypothesis that the Carina RV distribution 
is Gaussian can be discarded at 99% confidence level. The mean RV across the body of the galaxy 
varies from ^220 km s _1 at a distance of 7' (^200 pc) from the center to ^223 km s _1 at 13' (~400 
pc, 6cr level) and flattens out to a constant value of ~221 km s _1 at larger distances (600 pc, 4cr 
level) . Moreover and even more importantly, we found that in the Carina regions where the mean RV 
is smaller the dispersion is also smaller, and the RV distribution is more centrally peaked (i.e. the 
kurtosis attains larger values). The difference in mean RV is more than 4 km s" 1 (9<r level), when 
moving from E to W, and more than 3 km s _1 (~7c level), when moving from N to S. The RV gradient 
appears to be in the direction of the Carina proper motion. However, this parameter is affected by 
large uncertainties to reach a firm conclusion. There is evidence of a secondary maximum in RV across 
the Carina center (|£>|<200 pc). The occurrence of a secondary feature across the Carina center is 
also supported by the flat-topped radial distribution based on the photometric catalog. These findings 
are reminiscent of a substructure with transition properties, already detected in dwarf ellipticals, and 
call for confirmation by independent investigations. 

Subject headings: galaxies: individual (Carina) — galaxies: stellar content — galaxies: dwarf — Local 
Group — galaxies: kinematics and dynamics — techniques: radial velocities 
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1. INTRODUCTION 

Dwarf galaxies are fundamental laboratories to inves- 
tigate the influence of the environment on star forma- 
tion and on chemical evolution in stellar systems that are 
several order of magnitudes smaller than giant galaxies. 
Empirical evidence indicates that in the Local Group to- 
gether with the dwarf ellipticals (dE, M32-like) and the 
dwarf spirals (dS, M33-like) we are facing three differ- 
ent dwarf morphological types. The dwarf spheroidals 
(dSphs) show either single (Cetus, Monelli et al. 2010a; 
Tucana, Monelli et al. 2010b) or multiple star formation 
events (Carina, Bono et al. 2010), are spheroidal in mor- 
phology and devoid of neutral hydrogen (Bouchard et al. 
2005). The dwarf irregulars (dls) host a mix of old and 
young stellar populations, but they have recently experi- 
enced an intense SF episode (Sanna et al. 2009; Cole 
2009), moreover, they have irregular morphology and 
host a significant fraction of neutral hydrogen. The tran- 
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sition dwarfs have properties between the dSphs and the 
dls. The observational scenario concerning the dwarfs in 
the LG was recently enriched by the discovery of ultra- 
faint dSphs in the SDSS (Zucker et al. 2006; Tolstoy et 
al. 2009; Wyse 2010), but their properties still need to be 
investigated in detail. Although the morphological clas- 
sification appears robust, we still lack firm quantitative 
constraints concerning the evolution/transition between 
different dwarf morphological types. Kormendy (1985, 
1987) found a well defined dichotomy between ellipti- 
cals (Es) and spheroidals (Sphs, Kormendy et al. 2009). 
Early type galaxies are distributed along a sequence mov- 
ing from cD to dEs, while spheroidals are distributed 
along a sequence that overlaps with spirals and dwarf 
irregulars (see Fig. 1 in Kormendy et al. 2009). This 
empirical dichotomy seems to suggest that E and Sph 
galaxies arc stellar systems that underwent different for- 
mation and evolution processes. In particular, the Sphs 
might be either spirals or irregulars that lost their gas 
or transformed it into stars. However, there is no gen- 
eral consensus concerning the E-Sph dichotomy, since 
the correlation between the shape of the brightness pro- 
file and the galaxy luminosity is continuous when moving 
from Es to Sphs (Jerjen & Binggeli 1997; Gavazzi et al. 
2005; Ferrarese et al. 2006). 

Theoretical and empirical uncertainties affect our un- 
derstanding of several empirical correlations. More than 
twenty years ago it was suggested by Skillman et al. 
(1989) that the luminosity of dwarfs is correlated with 
their mean metallicity. However, it is not clear whether 
this correlation is linear over the entire metallicity range 
or shows a plateau in the metal-poor regime (Hclmi et 
al. 2006; Mateo et al. 2008; Kirby et al. 2008). More 
recently, it was suggested by Woo et al. (2008) that 
dEs/dSphs in the LG, at fixed total visual luminosity, 
show larger mass-to-light (M/L) ratios (see their Fig. 1). 
The same systems are, at fixed total stellar mass, more 
metal-rich than dls (see their Fig. 10 and Mateo 1998). 

Precise and homogeneous photometric and spectro- 
scopic data are required to address the above open issues 
concerning the formation and evolution of these systems. 
The advent of multi-object spectrographs at the 8-10m 
class telescopes provided the unique opportunity to col- 
lect sizable samples of low-, medium- and high- resolution 
spectra covering the entire body of nearby dwarf galaxies. 
By using hundreds of radial velocity (RV) measurements 
it has been suggested by Kleyna et al. (2004) that Sex- 
tans dSph host kinematic substructures, in particular, 
they found that the RV dispersion of the stellar popula- 
tion located across the center of this system is close to 
zero and increases outside the core. A different kinematic 
status has also been suggested by Battaglia et al. (2008) 
for the two distinct stellar populations in the Sculptor 
dSph. 

Moreover, evidence of a fall-off in the velocity disper- 
sion at large galactic distances was suggested in Sex- 
tans, Draco and Ursa Minor (Kleyna et al. 2003; Wilkin- 
son et al. 2004, 2006). However, Walker et al. (2009), 
using large and homogeneous samples of RV measure- 
ments based on high-resolution spectra collected with 
Mike Multi Fiber Spectrograph (MMFS) at Magellan 
(Walker et al. 2007, hereinafter W07), found that the RV 
profile in several dSphs (Carina, Draco, Leol, LeoII, Ursa 
Minor) are flat at large projected radii (R> 1 kpc) , while 



Fornax, Sculptor and Sextans show a gentle decline. 

The observational scenario is far from being settled, 
and indeed, Lokas et al. (2008) and Lokas (2009) using 
the same RV measurements provided by Walker et al. 
(2007), but a different algorithm to reject tidally stripped 
stars, found a significant decline in the RV dispersion 
profile of Carina, Fornax, Leol, Sculptor and Sextans. 
Furthermore, it was also suggested that kinematic status 
of Carina (Munoz et al. 2006) and Bootes (Belokurov et 
al. 2006) might be disturbed by the Milky Way. In this 
context, Carina plays a key role because it is relatively 
close (^=20.10 mag), shows at least two well separated 
SF episodes (t=l-6 and 12 Gyr) and it is metal-poor 
(<[Fe/H]> — 1.7 dex). 

Our group is involved in a long-term project on the evo- 
lutionary properties, variable stars, kinematic and abun- 
dances of stellar populations in the Carina dSph. This is 
the fourth paper of a series and it is focused on the ra- 
dial velocity distribution of candidate Carina stars. The 
structure of the paper is as follows. In Section 2, we 
discuss in detail the different spectroscopic data sets we 
collected for this experiment, together with the approach 
we adopted to reduce the data. Section 3 deals with ra- 
dial velocity measurements and with the observational 
strategy we adopted to validate radial velocities based on 
low-, medium- and high resolution spectra. In Section 4 
we lay out the selection criteria adopted to identify can- 
didate Carina stars and their radial velocity distribution. 
In Section 5, we investigate the velocity distribution, the 
velocity dispersion and the kurtosis of candidate Carina 
stars as a function of the projected radial distance. In 
this section, we also address the similarity between spec- 
troscopic and photometric radial distributions. Finally, 
in Section 6 we summarize the results of this investiga- 
tion and briefly outline possible future extensions of this 
photometric and spectroscopic experiment. 

2. SPECTROSCOPIC DATA SETS AND DATA REDUCTION 

To collect low-resolution spectra of old (Horizontal 
Branch [HB], V~20.75 mag), intermediate (Red Clump 
[RC], V~20.5 mag) and young (Blue Plume [BP], V~22 
mag) stellar tracers (see Fig. 1) we observed five fields 
with FORS2 at VLT (Appenzeller et al. 2008) in ser- 
vice mode 15 . The program was almost completed by 
one of us (M.N.) during a Visitor mode observations in 
December 2004. A significant fraction of FORS2 spec- 
tra were collected using the 1400V grism, which cov- 
ers Kp (A=4861.34 A) and Mgl Triplet (MgIT 16 ) lines, 
with a nominal resolution of 2100 (A c =5200 A). The slits 
used were 0"6 length, and typically 10" long, with a 
wavelength coverage «4560-5860 A. The exposure times 
range from 22 minutes (short) to more than one hour 
(long). 

During the December 2004 run, the 1028z grism was 
also used with a wavelength coverage «7730-9480 A. Its 
nominal resolution is 2560 (A c =8600 A) and covers the 
near- infrared (NIR) calcium triplet (CaT 17 ). Using the 
1400V grism we secured in 19.12 hrs of exposure time 743 
individual spectra of 356 stars with DIMM seeing rang- 

15 ESO programs 072.D-0671(A), P.I.: Bono; 078.B-0567(A), 
P.I.: Thcvcnin 

16 MgIT: A=5167.32, 5172.68, 5183.62 A 

17 CaT: A= 8498.02, 8542.09, 8662.14 A 
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ing from 0''35 to 1''27, while using the 1028z grism we 
secured in 1 hr of exposure time 113 spectra of 73 stars 
with seeing 0"38-0';59. Most of the stars with 1028z 
spectra also have 1400V spectra. We ended up with a 
total sample of 359 stars either with single or multiple 
FORS2 spectra. In particular, the LR data set collected 
with the grism 1400V includes 356 stars and among them 
eleven have a single spectrum, while the others a num- 
ber of spectra ranging from two to twelve. The LR data 
set collected with the grism 1028z includes 73 stars and 
among them 39 have single spectra, while the others a 
number of spectra ranging from two to three. In total the 
LR data set includes 359 stars of which eleven have a sin- 
gle spectrum and the others a number of spectra ranging 
from two to 14. Table 1 gives the log of the observa- 
tions of the spectra collected with FORS2. From left to 
right arc listed the date, the pointing, the coordinates, 
the grism, the exposure time and the seeing. 

These data were complemented by archival medium 
and high-resolution spectra of Red Giant (RG) and RC 
stars collected with GIRAFFE (Pasquini et al. 2002) at 
VLT 18 and covering the entire body of the galaxy (see 
Fig. 1). The GIRAFFE spectra (Koch et al. 2006) were 
collected using the grating LR8 with a nominal resolution 
of 6,500 centered on the NIR CaT. This data set (here- 
inafter GMR03) includes 10,394 individual fiber spectra 
of 1,070 stars and were secured with a total exposure time 
of 95.4 hrs with seeing condition of 0"4-2"2. Note that 
this data set was also adopted by Munoz et al. (2006) in 
their analysis of Carina radial velocity distribution. 

We also use the high-resolution spectra collected with 
the gratings: HR10 (R=19,800, 5339<A<5619 A), HR13 
(R=22,500, 6120<A<6405 A) and HR14A (R=17,740, 
6308<A<6701 A). This data set (hereinafter GHR) in- 
cludes 2,002 individual spectra of 98 stars, secured with 
a total exposure time of 24 hrs and seeing 0''57-3"20. 
The above data were complemented by a new set of GI- 
RAFFE spectra collected with LR8 and including 8,092 
individual spectra of 959 stars. This data set (hereinafter 
GMR08) was secured with 48.4 hrs of exposure time and 
seeing of 0"54-l"94. We ended up with a total sample 
of 1,931 stars either with single or multiple GIRAFFE 
spectra. 

The GMR (<GMR03+GMR08>, weighted mean) 
data set includes 1,887 stars and only one star has a single 
spectrum, while the others have multiple spectra ranging 
from two to 35. The HR data set includes 98 stars and 
all of them have multiple spectra ranging from seven to 
nine. Moreover, the GHR and the GMR03 spectra have 
38 objects in common, while 125 objects are in common 
between GHR plus GMR03 and GMR08 spectra. The 
LR resolution data set has 145 objects in common with 
the GIRAFFE (GHR plus GMR) data set. The objects 
in common were adopted to calibrate GMR and LR spec- 
tra. 

Data plotted in the top panels of Fig. 1 show the spatial 
coverage of the different data sets, while the bottom pan- 
els show the location of the targets in the V,B-I Color- 
Magnitude Diagram (CMD, Bono et al. 2010). The num- 
bers in parenthesis give the number of stars for which we 
estimated the radial velocity and the total number of 

18 ESO programs 074.B-0415(A), P.I.: Shctronc; 171.B-0520(B), 
P.I.: Gilmore; 180.B-0806(B), P.I.: Gilmore 



stars for which we collected at least one spectrum. 

Data plotted in the bottom panels of Fig. 1 show sev- 
eral interesting features: i) Stellar tracers - The FORS2 
spectra cover old, intermediate and young stellar trac- 
ers ( V<22, 359 RVs). ii) Statistics - Although the ex- 
posure time of the GIRAFFE spectra is a factor of 4-8 
longer than the FORS2 spectra, their limiting magnitude 
is V<20.6 (fibers vs slits). This means that this data set 
does not cover young (BP) and truly old (HB) tracers. 
The RGs are, indeed, a mix of old and intermediate-age 
stellar tracers. However, the sample of RV measurements 
based on GIRAFFE spectra (1,985 vs 359) is more than a 
factor of five larger than FORS2 (FoV ^25 vs ^5 arcmin 
squared). Hi) Multiplicity - A significant fraction of our 
targets have multiple spectra. This means that we can 
provide robust estimates of intrinsic errors, iv) - A sig- 
nificant overlap of low, medium and high-resolution spec- 
tra - We can constrain the occurrence of possible system- 
atic errors in RV measurements using spectra collected 
with different instruments and different gratings/grisms. 

We reduced our FORS2 data using standard IRAF 19 
tasks, using day time associated bias frames, flat field and 
calibration lamps. The wavelength calibration was car- 
ried out using the task identify, daytime lamp frames 
and a set of 8-12 calibration lines. We typically end 
up with an accuracy better than «0. 02-0. 04 A, <2.5 
km s _1 . However, we found systematic shifts up to 1 
pixel, corresponding to «30 km s _1 in radial velocity, 
for the two adopted grisms. To overcome the problem, 
we adopted the approach suggested by Kelson (2003). A 
model of the sky for each slit was created and subtracted 
from the reduced 2D spectrum. This model was also 
used to estimate and correct the systematic shifts from 
the daytime wavelength calibration lamps, using night 
sky lines, mainly OI 5577.34 A for the 1400V grism and 
a set of isolated skylines for the 1028z grism. 

Yet another potential source of systematics in RV esti- 
mates was the centroiding (see, e.g. Tolstoy et al. 2001). 
Thanks to the special care in the mask design, based 
on previous FORS2 preimaging of selected Carina fields, 
and to the very precise tracking of VLT, this effect was 
limited at the level of 5-10 km s _1 . To estimate the shift 
between the centroid of the star and the slit we adopted 
the slit image which is acquired just before the spec- 
trum. For each target the centroid of the slit is evaluated 
by collapsing the slit along its width and neglecting the 
pixels belonging to the star. The collapsed slit is then 
subtracted away from the pixels belonging to the star, 
whose centroid is then estimated with the task center 
in IRAF. It is worth mentioning that we also acquired slit 
images, for a couple of masks, soon after we collected the 
spectra. These slit images were adopted to confirm that 
our estimates of the centroids were minimally affected 
by tracking problems. The details of the pre-reduction 
and reduction of FORS2 spectra will be discussed in a 
forthcoming paper (Nonino et al. 2011, in preparation). 

The GIRAFFE raw data were retrieved from the ESO 
Archive. These data were reduced using IRAF. After 
the standard bias and flat correction, the spectra were 

19 IRAF is distributed by the National Optical Astronomy Ob- 
servatory, which is operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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extracted using the traces from flats. Wavelength cali- 
bration is based on daily calibration lamps: the formal 
solution rms was <0.02 A, corresponding to a system- 
atic <1 km s _1 in RV. After the fiber-spectra assign- 
ment, performed using the table associated with the raw 
data, all the spectra were visually inspected in order to 
remove bad spectra. A master sky for each exposure 
was created from sky fibers, after the spectra had been 
cross-correlated to remove shifts and scaled using the in- 
tensity of selected night sky lines in the wavelength range 
8250-8750 A, to better match the region with stellar ab- 
sorption lines. Subsequently the master skies were cross 
correlated with the targets spectra, scaled in intensity 
and subtracted. The cross correlation among the stacked 
skies showed that the relative accuracy in wavelength cal- 
ibration could give systematics of the order of 500 m s _1 
in RV. We finally coadded all the individual spectra of 
the same target using the task s combine, after correcting 
for the barycentric motion. 

The GIRAFFE spectra taken with the high resolution 
grisms were reduced following a similar approach. 

3. RADIAL VELOCITY MEASUREMENTS 

To measure the RV of GMR and GHR spectra we fit 
individual spectral lines in the coadded spectra. The 
RV of the LR spectra was measured on individual spec- 
tra and the final RV was estimated as a weighted mean 
among the multiple spectra, when available, of the same 
object. All spectra were normalized to the continuum 
and the fit to single or multiple lines performed using 
either a Gaussian or a Moffat function. The interactive 
code developed to perform the RV measurements will be 
described in a forthcoming paper (Fabrizio et al. 2011, 
in preparation). We selected by visual inspection sev- 
eral spectral lines for each data set, namely and Mgl 
triplet for FORS2/1400V spectra (see the spectra plot- 
ted in the panels d) of Fig. 2), the CaT for FORS2/1028z 
and GMR spectra (see the spectra plotted in the panels 
c) and b) of Fig. 2) and several strong Fel lines for GHR 
spectra (see the spectra plotted in the panels a) of Fig. 2). 
The quoted lines were typically fit with a Moffat func- 
tion with /3=2, to properly account for the contribution 
of the wings in the line fit. 

This preliminary estimate of the RV was validated 
by eye inspection and when judged satisfactory it was 
adopted to perform an automatic estimate of the RV 
using a large set of iron and heavy element lines. We 
adopted the line list for iron, a- and heavy-elements re- 
cently provided by Romaniello et al. (2008) and by Pedi- 
celli et al. (2010). The lines were selected according to 
the wavelength range and the spectral resolution of the 
different instruments and grisms. We ended up with a 
sub-sample of 30-70 lines for LR (grisms: 1028z, 1400V) 
spectra, with «35 lines for GMR spectra and with «90 
lines for the GHR spectra. For each spectrum we esti- 
mate the RVs as a weighted mean of the lines with the 
highest S/N ratio. This on average means 1-2 dozens 
of lines for LR and GMR spectra and approximately 40 
lines for the GHR spectra. 

The error in the RV based on individual lines was as- 
sumed equal to the sigma of the fitting function 20 . The 
RV of the entire spectrum was estimated as a weighted 

2(1 In particular, we assume eRv=c<J\/\ km s _1 



mean over the different fitted lines. The error in the ra- 
dial velocity of coadded spectra (GMR, GHR) is the error 
of the weighted mean. The RV of the LR spectra was 
estimated following the same approach for stars with a 
single spectrum and as a weighted mean among the dif- 
ferent mean measurements for the stars with multiple 
spectra. 

Data plotted in Fig. 3 show that the intrinsic error 
of the RV measurements based on GHR spectra is on 
average (biweight mean) 21 <1 km s _1 (panels a,b,c) and 
becomes of the order of 5 km s _1 for GMR03 spectra 
(panel d) . Note that the intrinsic error of GMR08 spectra 
is larger and —6 km s _1 . The difference is mainly due to 
the S/N ratio, since the latter sample is on average 0.5 
mag fainter and its total exposure time is almost a factor 
of two shorter. 

The intrinsic error of the RVs based on LR spectra 
is slightly larger and ranges from —11 km s _1 (grism 
1028z, panel f) to -10 km s" 1 ( grism 1400V, panel g). 
The difference between the FORS2 data sets is expected. 
The first two lines of the CaT, at the typical RV of Ca- 
rina stars, are contaminated by atmospheric OH emission 
lines. Therefore, the RV measurements based either on 
low- resolution spectra or on medium-resolution, low S/N 
ratio spectra might be affected by slightly larger uncer- 
tainties in case the sky lines are not precisely subtracted. 
The reader, interested in a more detailed discussion, is re- 
ferred to Nonino et al. (2007) and to Walker et al. (2007). 
Moreover, the 1028z spectra have a number of repeats 
that are approximately a factor of two smaller than the 
number of repeats of the 1400V spectra (see §2). 

The fraction of objects for which we measured the RV 
is quite high, and indeed the fraction ranges from —98% 
for GHR and GMR03 spectra, to «90% for LR (1400V) 
spectra and to more than «80% for GMR08 and LR 
(1028z) spectra. Typically the stars with no RV measure- 
ments have either low S/N or noisy spectra (see §2) or 
missing identification in the photometric catalog. Even- 
tually, the RV was measured in 2,165 out of the 2,323 
GIRAFFE spectra (-93%) and in 381 out of the 429 
(-89%) FORS2 spectra (see Fig. 3). In total the RV was 
measured in 1,812 stars using GIRAFFE spectra and in 
324 stars using FORS2 spectra (see Fig. 1). The final 
catalog includes RV measurements for 1,979 stars. 

To validate the approach adopted to measure the RVs, 
we decided to use real spectra. This approach has several 
indisputable advantages when compared with artificially 
generated spectra, since the different spectroscopic data 
sets we are dealing with partially overlap. We can con- 
strain the precision of RVs based on HR spectra, since 
for each object we have three spectra with similar spec- 
tral resolutions and covering three different wavelength 
regions (HR10, HR13, HR14A, see §2). The panels a) 
and b) of Fig. 4 show the internal comparison of our RV 
measurements as a function of the V-magnitude. The 
difference and the RV dispersion of the entire sample and 
of the candidate Carina stars (red dots, 180<RV<260 

21 Following the referee's suggestion, we adopted the biweight lo- 
cation estimator (Andrews et al. 1972), since it is a robust indicator 
insensitive to outliers in both Gaussian and non-Gaussian distri- 
butions. We adopted the definition given by (Beers et al. 1990), 
which includes data up to four standard deviations from the cen- 
tral location. This method is based on an iterative solution, but 
the process only requires a few steps to convergence. 
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km s _1 , see §4) are minimal. 

We can also constrain the precision of the GMR03 
spectra, since this subsample has more than three dozens 
of stars in common with the GHR spectra. The panel 
c) shows the difference between the RVs based on 
the <GHR> and those based on the GMR03 spectra 
(weighted mean). The mean (biweight) of the difference 
is 0.1 km s _1 for the candidate Carina stars (26) and 
smaller than 1 km s _1 for the entire sample (35). The 
GMR03 RVs were corrected for the difference. 

The panel d) shows the difference between the RVs 
based on <GHR> plus GMR03 (weighted mean) spec- 
tra and RVs based on GMR08 spectra. We have more 
than 120 stars in common and all of them are candidate 
Carina stars. The difference is smaller than -3 km s _1 
and the cr«8 km s _1 . The difference with the RVs based 
on GMR03 spectra is mainly due to the fact that the 
GMR08 spectra have a lower S/N ratio and are on av- 
erage fainter. The GMR08 RVs were corrected for the 
difference. 

The panel a) of Fig. 5 shows the difference between the 
RVs based on the GIRAFFE (<GHR+GMR>, weighted 
mean) spectra and on the FORS2- 1400V spectra (inter- 
nal comparison). The mean difference (biweight) over 
the entire sample of stars in common (140/144) is -3.6 
km s _1 and smaller (-2.7 km s _1 ) for the candidate Ca- 
rina stars (122/125), while the dispersion is of the or- 
der of 12 km s" 1 . The RVs based on FORS2-1400V 
were accordingly corrected. The panel b) shows the dif- 
ference between the RVs based on GIRAFFE and on 
FORS2 (<GHR+GMR> + 1400V, weighted mean) and 
RVs based on FORS2-1028z spectra. The difference in 
RV and dispersion are similar to the RVs based on the 
1400V grism and they were also corrected. Data plotted 
in panel c) show the difference between the RVs based on 
the GIRAFFE (<GHR+GMR>, weighted mean) and on 
the FORS2 (<1028z+1400V>, weighted mean) spectra 
after the corrections have been applied. The two samples 
have more than 140 stars in common and the difference 
is ~1 km s _1 for the entire sample and vanishing for can- 
didate Carina stars. 

The Fig. 6 shows the comparison between our RVs 
based on GIRAFFE spectra -<GHR+GMR>- with dif- 
ferent RV measurements available in the literature (ex- 
ternal comparison). The panel a) shows the comparison 
with the RVs based on echelle (+2D at LC02.5m, Mateo 
et al. 1993, crosses) and on multi-fiber spectra (HYDRA 
at CTIO 4m Blanco, Majewski et al. 2005, filled circles). 
The difference with the Mateo sample is quite small both 
in the mean (biweight) and in the velocity dispersion. 
The difference with the Majewski sample is larger (/i= 
3.6, cr=9.4 km s" 1 ), and it is caused by the lower S/N 
ratio of the spectra they adopted to estimate the RVs. 

The panel b) shows the comparison between our RVs 
and RVs based either on high-resolution spectra col- 
lected with MIKE at Magellan (crosses) or on medium- 
resolution spectra collected with GIRAFFE at VLT 
(filled circles) provided by Munoz et al. (2006). The 
MIKE spectra were collected in slit mode with a reso- 
lution of ~19,000 (red arm) and the RV measurements 
are based on the NIR calcium triplet. The GIRAFFE 
spectra adopted by Munoz et al. (2006) are the same 
GMR03 spectra we also included in this investigation. 
The difference with RVs based on MIKE spectra is min- 



imal (<1 km s _1 ), but the sample of stars in common is 
limited (14). The difference with the RVs based on GI- 
RAFFE spectra is larger ~3 km s" 1 if we account for the 
entire sample of stars in common (850/889 stars, (7=5.7 
km s _1 ) and for candidate Carina stars (353/373 stars, 
cr=6.7 km s _1 ). The reasons for this difference are not 
clear, apart from the fact that we adopted a different 
approach to prereduce the spectra and to measure the 
RVs. 

The panel c) shows the comparison with the RVs pro- 
vided by W0 7. The RV measurements provided by W0 7 
are based on high-resolution (20,000-25,000) spectra col- 
lected with MMFS at Magellan. These spectra sample 
the region across the magnesium triplet (5140-5180A). 
We have more than 550 candidate Carina stars in com- 
mon and the data plotted in this panel show that the 
two samples of RV measurements agree, within the errors 
quite well. The difference is -2 km s" 1 (biweight) and 
the RV dispersion is smaller than 10 km s _1 . The same 
outcome applies to the RVs based on FORS2 spectra (see 
panel d). We have more than 80 candidate Carina stars 
in common and the difference is wl km s _1 , while the 
dispersion er^ll km s . 

We ended up with a data set including RV measure- 
ments for 1979 stars. For the objects with RVs based 
on both GIRAFFE and FORS2 spectra we performed a 
weighted mean, using as weights the inverse square of the 
radial velocities standard errors of the GMR, GHR and 
LR spectra. As a final test, we compared our entire RV 
data set with the RVs provided by W07. Data plotted in 
panel e) further support the agreement between the two 
different sets of RV measurements, and indeed the differ- 
ence ranges from less than 1 km s _1 for the entire sample 
(785 objects in common) to 2 km s^ 1 for candidate Ca- 
rina stars (574), while the dispersion is ~9 km s _1 . The 
estimates of the mean difference are based, once again, 
on the biweight and the number of neglected objects is 
smaller than 9%. The precision of both internal and ex- 
ternal validations, and in particular the good agreement 
with the RVs provided by W07, support the approach we 
adopted to measure the RVs. 

The referee asked us to comment the impact of binaries 
stars on the current RV measurements. The occurrence 
of binary stars among RG stars in dSph galaxies is a 
highly debated topic and their impact on the RV dis- 
persion ranges from a sizable (Qucloz et al. 1995) to a 
small error (Olszewski et al. 1996; Hargreaves et al. 1996) 
when compared with the statistical error. However, in a 
recent investigation Minor et al. (2010) found, by using 
a detailed statistical approach, that dSph galaxies with 
RV dispersions ranging from 4 to 10 km s _1 can be in- 
flated by no more than 20% due to the orbital motion of 
binary stars. It is worth noting that the current LR spec- 
tra cover a time interval of three years (2004-2007) and 
the number of repeats ranges from two to 14 (only eleven 
stars with a single spectrum). The GMR spectra cover 
a time interval of five years (2003-2008) and the num- 
ber of repeats ranges from two to 35 (only 1 star with 
a single spectrum). The GHR spectra were collected in 
2005 and the number of repeats ranges from seven to 
nine. The three different data sets have 145 objects in 
common. This means that a significant fraction of stars 
in our sample have spectra collected on a time interval of 
several years. Therefore, the current RV measurements 
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are less prone to significant changes caused by binary 
stars. Moreover, the conclusions of the current investi- 
gation are minimally affected by a possible uncertainty 
of the order of 20% in the RV dispersion. 

4. CARINA RADIAL VELOCITY DISTRIBUTION 

The top panel of Fig. 7 shows the RV distribution of the 
entire (GIRAFFE+FORS2) data set. The well defined 
primary peak located at RV~220 km s" 1 includes candi- 
date Carina stars. The current radial velocity distribu- 
tion is soundly supported by the radial velocity distribu- 
tion recently provided by W07, but based on a smaller 
number of stars (see the middle panel of Fig. 7). This 
evidence and the minimal difference in the RV of the 
stars in common allowed us to merge the two RV cata- 
logs. For the objects in common weighted mean RVs were 
computed using as weights the inverse square standard 
errors in the individual RV measurements. We ended up 
with a sample of 2,629 stars (see the bottom panel of 
Fig. 7). 

To further constrain the radial velocity distribution of 
the candidate Carina stars we ran a Gaussian kernel on 
each star with a a equal to the RV uncertainty. The solid 
red line plotted in the top panel of Fig. 8 was computed 
by summing the individual Gaussians over the entire data 
set. We estimated the mean and the a by fitting the 
smoothed RV distribution with a Gaussian (dashed black 
line). We found that the peak (220.4±0.1 km s" 1 ) in 
the RV distribution agrees quite well, within the errors, 
with similar estimates available in the literature (Koch 
et al. 2006; Munoz ct al. 2006; Walker et al. 2007, 2009). 
The same outcome applies to the RV dispersion. Note 
that current estimate is based on a homogeneous sample 
of candidate Carina stars (180<RV<260 km s _1 , ^4er, 
1208 stars) that is ^55% larger than any previous sample 
of RVs. 

The middle panel of Fig. 8 shows the RV distribution 
of candidate Carina stars, but based on the RV measure- 
ments (780 candidate Carina stars) provided by W07. 
The two data sets provide, within the errors, very simi- 
lar kinematic properties. Note that by merging the two 
data sets of RV measurements we ended up, according to 
the quoted selection criterion, with a sample of 1378 can- 
didate Carina stars. This sample of Carina RVs is «75% 
larger than any previous RV sample. To overcome subtle 
uncertainties in the estimate of the RV due to the possi- 
ble presence of outliers in the final data set, we estimated 
the biweight mean of this sample and we found 220.9±0.1 
km s _1 . This estimate is based on 1369 stars and this 
is the sample that we will adopt in the following to esti- 
mate the RV distribution (red line in the bottom panel 
of Fig. 8). The mean based on the Gaussian fit (dashed 
line) agrees quite well with the biweight mean. However, 
the tails of the smoothed RV distribution appear shal- 
lower than the tails of the fitting Gaussian. There is also 
evidence that the low-velocity tail might be shallower 
than the high-velocity one. To provide a more quantita- 
tive estimate, we estimated the \ 2 of the two curves, and 
we found that the hypothesis that the smoothed distri- 
bution is Gaussian can be discarded at 99% confidence 
level. 

The referee has expressed concerns about the possible 
dependence of the Carina RV distribution on the different 
samples of RV measurements. To constrain this effect we 



estimated the RV distribution of Carina stars using stars 
brighter than V=20.25 mag (651). We fit the smoothed 
distribution with a Gaussian and we found a mean RV 
/i=221.3±0.1 km s" 1 and a a=8.7±0.1 km s" 1 . We 
also performed a new test including only stars brighter 
than V=20.50 mag (1050) and we found ^=220.9±0.1 
km s _1 and a <r=8.9±0.1 km s . The quoted values 
agree quite well with the estimates based on the entire 
sample. The two RV distributions based on smaller sam- 
ples of brighter objects also show asymmetric wings when 
compared with a Gaussian distribution. 

To further constrain the plausibility of the selection 
criterion to pin point candidate Carina stars, we plot- 
ted the kinematic candidates in the V,B-I (left panel of 
Fig. 9) and in the V, V-I (right panel of Fig. 9) CMDs. 
The symbols of the different spectroscopic data sets are 
the same as in Fig. 1. Data plotted in this figure show 
quite clearly the accuracy of the kinematic selection, and 
indeed the bulk of the candidate Carina stars are located 
along the expected evolutionary sequences (RG, RC, HB, 
blue plume) . The new sample also provides a clear sepa- 
ration between RG and Asymptotic Giant Branch (AGB, 
V~20, 5-/^1.6-1.9 mag) stars. Moreover and even 
more importantly, data plotted in Fig. 9 soundly support 
the results recently provided by Bono et al. (2010) con- 
cerning the metallicity distribution of Carina stars. The 
quoted authors found, using a color-color plane (U-V, 
B-I) to select candidate Carina stars, that the spread 
in color of RG stars is significantly smaller than sug- 
gested by spectroscopic measurements. The very narrow 
distribution in color of candidate Carina RGs is further 
supported by the kinematic selection. There are a few 
RGs that attains colors redder than typical Carina RGs. 
They might be either variables stars or misidentified ob- 
jects. 

5. DISCUSSION 

In order to constrain whether the occurrence of asym- 
metries in the radial velocity distribution is caused by 
the presence of substructures in Carina, we investigated 
the change in the RV distribution as a function of the 
projected radial distance (p). To avoid spurious fluc- 
tuations in the mean RV, we ranked all the candidate 
Carina stars as a function of p and estimated the mean 
(biweight) and the median using the first 200 objects. 
The p of each subsample was estimated as the mean over 
the individual distances of the same 200 stars. We es- 
timated the same quantities by moving of one object in 
the ranked list until we accounted for the most distant 
200 stars in our sample located inside the tidal radius 
(r t ~28.8±3.6', Mateo 1998). The error on the mean RV 
for individual bins is smaller than one tenth of km s . 
In order to provide robust constraints on the possible un- 
certainties introduced by the number of stars per bin and 
by the number of stepping stars we performed a series of 
MontcCarlo simulations. The estimated mean dispersion 
is plotted as cyan shaded area across the RV mean in the 
bottom panel of Fig. 10. Note that changes in the binning 
criteria can change the extent of the secondary features, 
but their relative positions are minimally affected. 

A glance at the data plotted in the bottom panel 
of Fig. 10 shows that the mean RV attains a well 
defined minimum -RV~219.5±0.4 km s" 1 - for p -7' 
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(—200 pc) 22 from the galaxy center and increases up to 
-222.6±0.4 km s" 1 for p -13' (-400 pc, 6a level) and 
a plateau of — 221.4±0.4 km s _1 at larger distances (600 
pc, 4<t level). This evidence is robust, since the mean 
(biweight, black line) and the median (blue line) attain 
very similar values and intrinsic errors. To further con- 
strain the change of the kinematic status across the body 
of the galaxy we also estimated the radial velocity disper- 
sion (crRy, solid line middle panel of Fig. 10). Note that 
to avoid subtle uncertainties caused by asymmetric ra- 
dial velocity distributions we also estimated for each bin 
the semi-interquartile range (SIQ, dashed-dottcd line in 
the middle panel of Fig. 10). Data plotted in the middle 
panel of Fig. 10 show that both the gev and the SIQ 
show very similar radial trends, but the former disper- 
sion parameter attains values that are «50% larger. The 
difference indicates that the velocity dispersion should be 
cautiously treated in dealing with asymmetric RV distri- 
butions (Strigari et al. 2010). We also found that the 
dispersion is larger in the innermost regions, decreases 
in coincidence of the minimum in the mean RV and after 
a mild increase attains an almost constant value in the 
outermost galaxy regions. Note that the cyan shaded 
areas across the RV dispersion and the SIQ display the 
intrinsic error on individual bin according to MonteCarlo 
simulations. There is a hint that the dispersion could in- 
crease in the outskirts, thus supporting the finding by 
W07. However, more data are required to constrain the 
change in the projected radial velocity dispersion across 
the tidal radius. 

Interestingly enough, the kurtosis plotted in the top 
panel of Fig. 10 shows a mirror trend compared with the 
projected radial velocity dispersion and with the mean 
RV. This evidence is suggesting that in the Carina re- 
gions where the mean RV is smaller, the dispersion is 
also smaller and the radial velocity distribution more 
centrally peaked. The opposite trend characterizes the 
regions in which the mean RV is larger. The outermost 
regions show fiat trends in the quoted parameters. Evi- 
dence of kinematic and photometric substructures have 
already been found in several LG dwarfs (Monelli et al. 
2003; Kleyna et al. 2004; Walker et al. 2009). In partic- 
ular, Battaglia et al. (2008) found evidence of a radial 
velocity gradient in Sculptor, that they interpreted as an 
indication of internal rotation. The current data indi- 
cate that the RV distribution across the body of Carina 
shows a radial gradient and possible evidence of rotation. 
However, the RV distribution is far from being uniform, 
possibly suggestive of the occurrence of well defined spa- 
tial substructures. 

The referee has expressed concerns about the approach 
we adopted to estimate the weighted mean between the 
GIRAFFE and the FORS2 spectra. To clarify this key 
point, Fig. 11 shows the same data of Fig. 10, but the pa- 
rameters plotted in the left panels are only based on the 
GIRAFFE spectra, while those plotted in the right pan- 
els are based on the weighted mean between GIRAFFE 
and FORS2 spectra. The trends of the plotted param- 
eters as a function of the radial distance are the same. 

22 Projected distance on Carina were estimated assuming a 
true distance modulus of 20.10±0.12 mag (DalPOra et al. 2003; 
Pietrzyriski et al. 2009) together with the inclination angle and the 
axes ratio (Mateo 1998). 



They are only slightly more noisy due to the decrease in 
the sample size. 

To further constrain the spatial extent of possible kine- 
matic substructures we investigated the same parame- 
ters, but they were estimated as a function of right ascen- 
sion (a) and declination (8). Data plotted in panel c) of 
Fig. 12 show the marginal of the radial velocity along the 
right ascension axis. We found that the mean RV, when 
moving from E (—750 pc from the galaxy center) to W 
(—-1060 pc from the center), decreases from 221.9±0.4 
to 219.8±0.4 km s _1 . The difference is equal to -2.1 
km s _1 (4<T level) and increases by almost a factor of two 
(—4.5 km s _1 , 9a level) if we take into account the abso- 
lute maximum (— 223.3±0.4 km s _1 , —320 pc from the 
center) and the absolute minimum (— 218.8±0.4 km s _1 , 
— -280 pc from the center) of the mean RV. 

It is worth mentioning that the increase in the mean 
velocity both in the external regions and outside the 
galactic center appears to be in the same direction of 
the Carina proper motion (large blue arrow in panel d) 
of Fig. 12, Piatek et al. 2003; Walker et al. 2009). How- 
ever, the uncertainty affecting current estimates of the 
Carina proper motion is quite large (see small blue ar- 
rows). Note that the increase in the mean RV cannot 
be in the direction of the Galactic center (heavy dashed 
line). Moreover and even more importantly, we found 
that in a region of 200 pc across the galactic center the 
RV shows a secondary maximum that is strongly corre- 
lated with the velocity dispersion (panel b) and anticor- 
related with the kurtosis (panel a). 

The change in the mean RV along the declination axis 
is similar to the change in the right ascension axis. The 
mean RV decreases, when moving from N to S (panel e), 
from 221.9±0.4 (-630 pc from the center) to 219.8±0.4 
km s _1 (—-650 pc from the center). The difference is at 
least at 4a level, while the difference between maximum 
and minimum (222.3±0.4 km s _1 , -238 pc; 219.1±0.4 
km s _1 , —41 pc) is more than 3 km s _1 (— 7a level). 
The absolute minimum is located close to the center in 
the IV quadrant and associated to the absolute mini- 
mum in RA. The radial trend outside this region shows 
the quoted secondary maximum and a constant smooth 
trend at larger distances from the center. The dispersion 
(panel f) and the kurtosis (panel g) are once again corre- 
lated/anticorrelated with the behavior of the mean RV. 
It is noteworthy, that the quoted findings are minimally 
affected by the FORS2, GIRAFFE and W07 data sets 
adopted to estimate the RVs of candidate Carina stars. 

To explain the occurrence of a kincmatically cold popu- 
lation close to the center of Sextans, Kleyna et al. (2004) 
suggested, using also photometric radial distributions, 
that it could be the aftermath of a merging with a glob- 
ular cluster. To further constrain the nature of the above 
kinematic peculiarities we investigated the Carina radial 
distributions, using the same approach adopted by Mon- 
elli et al. (2003) and the recent photometric catalog by 
Bono et al. (2010). The isocontour levels (panel b) and 
the marginal (panel a) along the right ascension axis plot- 
ted in Fig. 13 show that the radial distribution is flat- 
topped across the Carina center. Note that the bright 
field star HD48652 ( V=9.14 mag) minimally affects this 
feature, since it is located at 2'77 from the Carina center. 
On the other hand, the radial distribution is symmetric 
and with a well defined peak along the declination axis. 
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We drew three vertical lines across the center and the two 
secondary maxima. Interestingly enough, the same lines 
include the secondary maximum detected along the right 
ascension axis in the mean RV. These empirical evidence 
indicates that the regions across the galaxy center might 
form a substructure probably reminiscent of a transition 
between a bulge-like and/or a disk-like substructure. 

6. SUMMARY AND CONCLUSIONS 

We presented new radial velocity measurements of 
Carina dSph stars based on low-resolution spectra 
(R/-2500) collected with the FORS2 multi-object slit 
spectrograph at the VLT. The key advantage of the cur- 
rent sample is that for the first time we collected spectra 
of old (HB, V-20.75 mag), intermediate (RC, K~20.5 
mag) and young (BP, K~22 mag) stellar tracers. The 
bulk of the spectra were collected using the 1400V grism 
covering the Hp and the Mgl triplet. For a fraction (20%) 
of these targets we also collected spectra using the 1028z 
grism centered on the NIR CaT. We secured 856 spectra 
of 359 stars with ^20 hours of exposure time. 

The above data were complemented by medium 
(R~6500) and high-resolution (R~18000-23000) spectra 
collected with the GIRAFFE multi-object fiber spectro- 
graph at the VLT. The targets of this spectroscopic data 
set are RG and RC stars (V<20.6 mag). The medium 
resolution spectra are centered on the NIR CaT, while 
the high-resolution ones cover a broad wavelength region 
(5300<A<6700 A). In total this data set includes 20,488 
spectra of 1931 stars collected with ~168 hours of expo- 
sure time. 

The radial velocity distribution based on the above 
spectra agrees quite well with similar data available in 
the literature. By assuming candidate Carina stars the 
objects with RVs ranging from 180 to 260 km s _1 (4<r) 
we ended up with a sample of ~1,210 stars. By fitting 
these objects with a Gaussian we found a mean RV - 
<RV>=220.4±0.1 km s _1 - and a velocity dispersion - 
<r=11.7±0.1 km s _1 - that agree, within the errors, with 
similar measurements provided by Walker et al. (2007). 

The above kinematic selection was firmly supported by 
the position of the candidate Carina stars along the ex- 
pected evolutionary sequences (RG, AGB, RC, HB, Blue 
Plume) of the CMD. The narrow distribution in color 
of the candidate Carina RG stars confirms the recent 
findings by Bono et al. (2010), based on photometric se- 
lection criteria, that the spread in metallicity of Carina 
stars is smaller than suggested by spectroscopic measure- 
ments (Koch et al. 2006). 

The current catalog of Carina RV measurements was 
complemented by similar measurements recently pro- 
vided by Walker et al. (2007). The stars in common were 
averaged and we ended up with a sample of ~ 1,3 70 RV 
measurements of candidate Carina stars that is «75% 
larger than any previous Carina RV sample. We found 
that the hypothesis that the Carina RV distribution is 
Gaussian can be discarded at the 99% confidence level. 
A more detailed investigation of the RV across the body 
of the galaxy indicates that the mean RV changes from 
^220 km s _1 at a distance of 7' (~200 pc) from the cen- 
ter to ^223 km s _1 at 13' (6<r level) and attains an al- 
most constant value of ~221-222 km s _1 at larger dis- 
tances (3-5cr level). Moreover, we found that in the Ca- 
rina regions where the mean RV is smaller the dispersion 



also attains smaller values and the RV distribution is 
more centrally peaked (larger kurtosis). 

The mean RV at large distances from the galactic cen- 
ter decreases, when moving from E to W, by more than 
2 km s^ 1 (4cr level). The difference increases by almost 
a factor of two (^4 km s _1 , 9a level) if we account for 
the difference between the absolute maximum and the 
absolute minimum of the mean RV. The change in the 
mean RV, when moving from N to S, is similar. The dif- 
ference increases from more than 2 km s _1 (4a level) in 
the outermost regions to more than 3 km s _1 (la level) 
across the galaxy center. The RV gradient appears to be 
in the direction of the Carina proper motion. However, 
this parameter is still affected by large uncertainties to 
reach a firm conclusion. 

Furthermore, there is also evidence of a secondary max- 
imum in RV across the Carina center (|D|<200 pc). The 
occurrence of a secondary feature across the Carina cen- 
ter is also supported by the flat-topped radial distribu- 
tion (photometric catalog) along the right ascension axis. 
In particular, the two secondary maxima cover the same 
galaxy regions of the secondary maximum in the RV ra- 
dial variation. These results, once independently con- 
firmed, are probably reminiscent of a substructure with 
transition properties. 

Evidence of substructures in dwarf galaxies dates back 
to Jerjen et al. (2000) who detected a spiral structure 
in the Virgo dE IC3328 using deep, ground-based R- 
band images. This discovery was soundly confirmed by 
Liskcr et al. (2007) who found that a significant fraction 
of bright (Mb < 16) early-type dwarfs in Virgo are char- 
acterized by disk- like features (spiral arms, bars). Fur- 
thermore, they also found that the properties of dE in 
clusters are also correlated with environmental density. 
This evidence supports numerical simulations of galaxy 
harassment indicating that late-type galaxies undergo 
a significant transformation when accreted in a cluster 
(Mastropietro et al. 2005; Kormendy et al. 2009). It was 
also suggested that dEs with disk-like features might be 
the low-luminosity tail of normal disk galaxies (Lisker 
et al. 2006). This working hypothesis was recently sup- 
ported by new and accurate kincmatical data of IC3328 
suggesting that the observed velocity dispersion is the 
aftermath of two distinct substructures, namely a thin 
stellar disk and a dynamically hot component. However, 
no robust conclusion could be reached due to limited ra- 
dial coverage (Lisker & Fuchs 2009). 

Irregular kinematic properties have also been found 
in the recently discovered ultra-faint dwarf galaxy - 
Willman I- by Willman et al. (2010) using both photo- 
metric and spectroscopic data. They found that galaxy 
stars located in innermost regions show a radial velocity 
offset of 8 km s _1 when compared with the outermost 
ones. Moreover, they also found initial hints of asymme- 
tries in the radial velocity distribution, but their spec- 
troscopic sample is too small to reach firm conclusions. 

If the complex kinematic properties of Carina will be 
supported by future more accurate and deep data sets 
then the occurrence of disk-like features might be consid- 
ered a wide spread property of gas-poor dwarf galaxies. 
A detailed comparison between observed and predicted 
radial velocity distribution and radial velocity dispersion 
profile is required to constrain the nature of Carina kine- 
matics. 
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Fig. 1. — Top Left: Spatial distribution of different spectroscopic data sets. The blue pluses and the red triangles mark targets with 
medium (LR03,LR08) and high-resolution (HR10,HRf 3,HR14A) GIRAFFE spectra. The background (grey dots) is a composite Carina 
reference image based on randomly selected subsamples of stars from the Bono et al. (2010) photometric catalog. The red dashed ellipse 
shows the inclination (i=65°) and the radial extent (tidal radius, rt=28.8±3.6' Mateo 1998). Top Right: Same as the top left, but for 
targets with low-resolution FORS2 spectra (green diamonds). The yellow asterisks mark the positions of the targets with high-resolution 
spectra collected with MMFS at Magellan by Walker et al. (2007). Bottom Left: V ,B-I CMD of spectroscopic targets. The colored 
symbols are the same as in the top left panel and show the position in the CMD of the targets with GIRAFFE spectra. The numbers in 
parentheses are referred to the number of stars with measured radial velocities and to the total number of stars with at least one observed 
spectrum. Bottom Right: Same as the bottom left, but for targets with FORS2 and MMFS spectra. 
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Fig. 2. — Selected spectra collected with different spectrographs, a): High-resolution spectra of three RGs, with apparent V magnitudes 
ranging from ~17.6 to ~19.3, collected with GIRAFFE at VLT. The data plotted in the three panels show a limited wavelength range 
(5428-5461 A) around four strong Fel lines (vertical dashed lines). The individual radial velocities (RV) are also labeled. The numbers in 
parentheses arc referred to the IDs of individual stars in the photometric catalog, b): Same as panel a), but for medium- resolution spectra 
collected with GIRAFFE at VLT. These spectra are centered on the NIR calcium triplet (CaT). The top and the middle panel show the 
spectrum of two RGs, while the bottom panel the spectrum of a RC star (y~20.5 mag), c): Same as panel a), but for low-resolution 
spectra collected with FORS2 at VLT. These spectra are centered on the calcium triplet (CaT, grism: 1028z). The top panel shows the 
spectrum of a RG star (V~17.8 mag), the middle panel the spectrum of a RC star (V~20.5 mag) and the bottom panel the spectrum 
of a HB star (V~20.8 mag), d): Same as the panel c), but for low-resolution spectra centered on the magnesium triplet (MgIT, grism: 
1400V). The top panel shows the spectrum of a RG star ( V~18.0 mag), the middle panel the spectrum of a RC star ( V~20.6 mag) and 
the bottom panel the spectrum of a HB star (l/~20.8 mag). 
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Fig. 3. — Intrinsic errors in RV measurements of different spectroscopic data sets as a function of visual magnitude. From top to bottom 
RV errors based on GHR (panels a,b,c), on GMR (panels d,e) and on LR (panels f,g) spectra. In each panel the labels give the wavelength 
range, the resolution of the spectra and the biweight mean of the intrinsic error. The number of RV measurements and the total number 
of spectroscopic targets are also labeled in the top right corner. 
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Fig. 4. — Comparison of RV measurements based on medium- (GMR) and high-resolution (GHR) spectra collected with GIRAFFE at 
VLT as a function of the visual magnitude. The panels a) and b) show the difference in RV (ARV=RVio-RVi3 i i4) among the GHR spectra. 
The biweight mean (fi), the standard deviation for the entire sample (black labels) and for the candidate Carina stars (180<RV<260 km s _1 , 
red labels) are also labeled. The numbers in square parentheses show the number of objects in common between the two samples before 
and after the biweight mean. The panels c) and d): Same as top panels, but the difference is between the weighted mean of the radial 
velocity based on the entire sample of GHR spectra and on the GMR spectra. 
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Fig. 5. — Same as Fig. 4, but the difference is between the weighted mean of the RV measurements based on all the GIRAFFE 
(GHR+GMR) spectra and those based on the FORS2-1400V low-resolution (LR) spectra ( ARV =RV gh R+GM fl-RV FORS— 1400V > panel 
a). The panel b) Same as panel a), but the difference is between RVs based on GIRAFFE plus FORS2-1400V spectra (weighted mean) and 
those based on the FORS2-1028z spectra. The panel c) - Same as panel a), but the difference is between the entire sample of GIRAFFE 
and FORS2 RV measurements. The labels and the numbers in parentheses have the same meaning of Fig. 4. 
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Fig. 6. — Comparison among RV measurements based on different spectroscopic data sets. The panel a) shows the difference between our 
GIRAFFE (GHR+GMR) RVs and those provided by Mateo et al. (1993) and by Majewski et al. (2005). The panel b) shows the difference 
between our GIRAFFE (GHR+GMR) RVs and those provided by Munoz et al. (2006) using two different data sets. The panel c) and d) 
show the difference between our GtRAFFE (GHR+GMR) and FORS2 RVs with those provided by W07 (Walker et al. 2007) using spectra 
collected with MMFS at Magellan. The panel e) shows the difference between our entire sample (GIRAFFE+FORS2) of RVs and those 
provided by W07. The biweight mean (fi), the standard deviation for the entire sample (black labels) and for the candidate Carina stars 
(180<RV<260 km s , red labels) are also labeled. The numbers in square parentheses show the number of objects in common between 
the two samples before and after the biweight mean. 
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Fig. 7. — Top: Radial velocity distribution of the entire sample normalized to the maximum as a function of the radial velocity. Middle: 
Same as the top, but based on RV measurements provided by W07 (Walker et al. 2007). Bottom: Same as the top, but based on both our 
and W07 RV measurements. Note that for the stars in common in the two data sets we computed a weighted mean. See text for details. 
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Fig. 8. — Logarithmic radial velocity distribution for the candidate Carina stars (180<RV<260 km s -1 , ~4<x). The red solid line shows 
the smoothed radial velocity distribution estimated running a Gaussian kernel on individual RV measurements. The mean (/i) and the <r 
of the Gaussian fit (dashed black line) are labeled together with the total number of candidate Carina stars. Middle: Same as the top, 
but based on RV measurements provided by W07 (Walker et al. 2007). Bottom: Same as the top, but based on both our and W07 RV 
measurements. 
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Fig. 9. — V,B-I(\eft) and V, V--T(right) CMDs of candidate Carina stars selected on the basis of the radial velocity (180<RV<260 
km s — , ~4cr). The symbols of the different spectroscopic data sets are the same as in Fig. 1. 
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Fig. 10. — Bottom: Mean (biweight, black line) and median (blue line) projected mean radial velocity (km s -1 ) as a function of radial 
distance in arcminutes (bottom axis) and in parsecs (top axis). The individual values were estimated as a running average over sub-samples 
of 200 stars (see text for more details). The cyan and the purple shaded areas across the mean and the median curves display a generous 
estimate of the intrinsic error on individual bins according to MonteCarlo simulations. Middle: Same as the bottom, but for the projected 
radial velocity dispersion (o\rvj black line) and the semi-interquartile range (SIQ, blue line). The cyan shaded areas across the curves 
display the intrinsic error. Top: Same as the middle, but for the kurtosis. 




Fig. 11. — Left — Same as Fig. 10, but the estimates of the different parameters are based only on GIRAFFE (GMR+GHR) spectra. 
Right — Same as the left panels, but the estimates are based on the weighted mean between GIRAFFE and FORS2 spectra. 
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Fig. 12. — The panels a),b),c)- Same as Fig. 10, but as a function of right ascension (a, bottom axis) and distance (pc, top axis). The 
panel d)- Sky distribution of the different spectroscopic data sets. Symbols and colors are the same as in Fig. 1. The long blue arrow shows 
the Carina proper motion according to Piatek et al. (2003) and to Walker et al. (2009), while the short ones display current uncertainties. 
The dashed black line shows the direction of the Galactic center according to an observer located in the center of Carina. The dotted and 
the dashed-dotted thin lines display the secondary features identified in Fig. 10 and in Fig. 13. The panels e),f),g)- Same as Fig. 10, but 
as a function of declination (<5, left axis) and distance (pc, right axis). 
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Fig. 13. — The panel a)- Normalized radial profile along the right ascension axis of the Carina photometric catalog (Bono et al. 2010). 
The panel b)— Sky distribution of the photometric catalog. Red contours display iso-density levels ranging from 5 to 95% with steps of 
~11%. The two thin dashed-dotted lines mark the two secondary peaks, while the thin dotted line marks the center of the galaxy. The 
blue star shows the position of the bright field star HD48652 ( F=9.14 mag). The panel c)— Same as the panel a), but along the declination 
axis. 
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Right Ascension (J2000), the units arc hours, minutes and seconds. 
^Declination (.T2000), the units arc degrees, arcminutcs and arcseconds. 
C Exposurc time (seconds). 

^Initial and final DIMM seeing (arcseconds) during the exposure time. 

C The pointings f3 and f3r arc the same, but the latter was rotated by 180 degrees (Nonino ct al. 
2011, in preparation). 



